Dynamics of anisotropic spin dimer system in strong magnetic field 
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Recently measured high-field ESR spectra of the spin dimer material TlCuCla are described within 
the framework of an effective field theory. A good agreement between the theory and experiment is 
achieved, for all geometries and in a wide field range, under the assumption of a weak anisotropy of 
the interdimer as well as intradimer exchange interaction. 
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I. INTRODUCTION 

Gapped spin systems in high magnetic field have at- 
tracted much attention recently, both from the theoret- 
ical and experimental side. In absence of the field the 
system is supposed to have a singlet ground state and 
a finite gap A to the lowest excitation, typically in the 
triplet sector. When the field is increased beyond the 
critical value He necessary to close the gap, the ground 
state acquires a finite magnetization, and a number of 
new phenomena can appear, including critical phases, 
field-induced ordering, magnetization plateaux, etc. 

The system behavior aX H > He depends strongly on 
the symmetry properties and dimensionality. Generally, 
if there is no axial symmetry with respect to the field di- 
rection, the high-field phase always exhibits a tranverse 
staggered long-range order (LRO) perpendicular to the 
applied field, independently on the system dimensional- 
ity, and is characterized by a finite spectral gap. If the 
axial symmetry is present, it cannot be spontaneously 
broken in a one-dimensional (Id) system. The high-field 
phase is in this case characterized by quasi-LRO (power- 
law correlations), and its lowest excitations are deter- 
mined by the spinon continuum. 

In the 3d case, U(l) symmetry gets spontaneously bro- 
ken and the high-field phase is ordered but possesses a 
gapless (Goldstone) mode. If one views the process of for- 
mation of the high-field phase as accumulation of hard- 
core bosonic particles (magnons) in the ground state, the 
ordering transition aX H — He can be interpreted as the 
Bose-Einstein condensation (BEG) of magnons. The idea 
of field-induced BEG was discussed theoretically many 
timesji The best available realisation of such a transi- 
tion was observe d^i^'^ in TlGuGla, which can be viewed 
as a system of coupled S — \ dimers. The temperature 
dependence of the uniform magnetization was found^ to 
agree qualitatively with the BEG theory predictions. 

High-frequency magnetic resonance measurements^ 
have demonstrated directly the field dependence of the 
energy gap, but were limited to the fields H < He and 



a few number of microwave frequencies. The response 
in the high-field ordered phase of TlGuGla was measured 
in the inelastic neutron scattering (INS) experiments of 
Riiegg et al.;:'''4 the behavior of the lowest triplet gaps 
as functions of field was successfully described within 
the bond-operator mean- field theory.^ In those INS mea- 
surements only two triplet modes were observed above 
He, and the gap of the third mode was concluded to 
be zero within the experimental resolution; however, the 
low-energy range could not be studied in this experiment 
because of the strong field-induced magnetic Bragg con- 
tamination below 0.25 meV (60 GHz). 

Very recently, high-field electron spin resonance ex- 
periments on TlCuCla in a wide range of fields up to 
90 kOe were conducted, which revealed a reopening of 
the gap above He in the low-energy range inaccessible by 
means of INS. A natural explanation would be the exis- 
tence of some anisotropic interactions explicitly breaking 
the U(l) symmetry. The aim of the present work is to 
show that the available data can indeed be described on a 
quantitative level, assuming presence of a weak exchange 
anisotropy in both intra- and interdimer interactions. 



II. EXPERIMENTAL SUMMARY 

The crystals of TlGuGla have monoclinic symmetry, 
with crystallographic axes a and c forming an angle of 
96.32°, b being the twofold axis. The sample growth is 
described in detail in Ref. 0- ESR spectra were taken at 
temperature 1.5 K, in the field range from to 90 kOe, 
using a set of home-made microwave spectrometers with 
transmission type cavities and a superconducting mag- 
net. Single crystals with the volume of 20 50mm^ were 
used. During the experiments crystals were mounted in 
the following orientations with respect to the magnetic 
field: -ff||[010] (i.e., parallel to the b axis), H _L (102) 
and i/||[201]. The [201] direction forms an angle of 51.7° 
with the a axis, see Fig. ^ The detailed description of 
the ESR spectra is reported in Ref. and here we will 
just summarize the results. 
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FIG. 1: Schematic picture of the anisotropy axes in TlCuCla. 

The ESR signal corresponding to the transitions from 
the ground state to the lowest g = excited state is 
clearly visible both below and above the critical field. 
For H < He, such transitions would be forbidden for an 
ideal isotropic system. The main feature of this signal is 
reopening of the gap above the critical field, which can 
be directly observed for all three mutually perpendicu- 
lar field directions. Besides the ground state transitions, 
transitions between Zeeman-split components of the ther- 
mally activated triplet were observed below the critical 
field. The analysis of the field dependence of thermally 
activated transitions aX H < He has also suggested pres- 
ence of a finite zero-field splitting of the triplet >i2i 

III. THE EFFECTIVE MODEL 

The dynamics of a 3d coupled anisotropic -S* = ^ dimer 
system in a wide range of fields can be described within 
the effective field theorjsiiii^ which may be viewed as a 
continuum version of the bond boson approachiSiH It is 
based on introducing dimer coherent statesii 

|A, B) = {l-A'- B^Y/^\s) + Y,{A, + (I) 

i 

where \s) and j = {x,y,z) are the singlet and three 
triplet statesji^ and A, B are real vectors related to the 
magnetization M = {Si + S2) and sublattice magneti- 
zation L = {Si — S2) of the spin dimer: 

M = 2{AxB), L = 2A^/l- - B^. (2) 

In TlCuCla, there are two types of chains of spin dimers 
running along the crystallographic a axis. In the high- 
field phase the staggered order L alternates between the 
two different types of chains.** The magnetization M 
should be uniform, which leads us to the following ansatz: 

where R\^i, = Aei -I- /xe2 -I- ve^ is the radius- vector label- 
ing the dimers, and the vectors ei 2,3 are in the following 
way connected to the lattice vectors a, b, c: 

ei^a, e2 = (6-c)/2, e3 = (6 + c)/2. 

The topology of exchange paths in TlCuCla is known 
rather well.^'l^ We assume additionally a weak or- 
thorhombic anisotropy (the simplest type compatible 



with the crystal symmetry) in intra- as well as in inter- 
dimer exchangejiS< so that instead of one intradimer ex- 
change constant J one has a vector J = {Jx, Jy, Jz}, etc. 
Each dimer has an inversion center, which excludes the 
intradimer Dzyaloshinskii-Moriya (DM) interaction (but 
it remains possible between the dimers); for simplicity, 
we neglect the DM interaction in the present treatment. 
Passing to the continuum, one obtains the Lagrangian 

C = -2hri ■dt^-\ J2 E [(^^- • '^)^*]' (3) 

i—x,y,z j — 1.2,3 

- "^{m-^^p^ + miijf} + 2h-{<fXT])- V{ip, rj), 

i 

where rrii — fhi — fii and fhi = \\(-ijn\{Jj + Jn)', we de- 
note hi — gijfXBHj, where g is the gyromagnetic ten- 
sor and H is the magnetic field. Effective interdimer 
couplings Pi are in the following way connected with the 
microscopic couplings (the notation of Ref. ^^is used): 

/3 = «^(ioo) + '^(201) + ^'^(iH) ~ ^''(100) ~ 2J(\i 1). (4) 

We will assume that we are not too far above the criti- 
cal field, so that the magnitude of the triplet component 
is small, i.e., ip,ri <^ 1. Then, retaining only the fourth- 
order terms in the interaction 1^ in one obtains 

i ij 

lij = Pi - Aij, Kj = - ^ \f-iji\ai, (5) 

where ai is given by a different combination of couplings, 
" = '^(loo) + ^^(201) + 2^(111) + 2J(ioo) + 2J('^i). (6) 

The further derivation is similar to that given in Ref. 
[1^ Integrating out the field 77 results in the relation 

'Ht = QijFj, F = -Wt^p + {hx (f) 

dij ^ 6^jLpI ip^ipj 

Wij = — / ,lki ^ij ~ ~ • {() 

rrii rriiraj 
and yields the effective Lagrangian depending on only: 

i i—x.y,z J — 1,2,3 

- y^2^ihx^),dm-U2{<p)-U4{<p,dt<p), (8) 
^ — ' rrii 

i 

with the quadratic and quartic interaction given by 

U2 = m,^^^ - ^{h X (9) 
mi 

^-^ ^-^ mi ^-^ miTrij 

i ij * ij 
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We analyze the obtained field theory at the mean-field 
(zero-loop) level; this is formally justified for weak inter- 
dimer coupling, while in TlCuCla /?; and rhi are of the 
same order. However, a similar approximation is used 
in the bond-boson approach which successfully describes 
the INS data on the magnon dispersion in this materiali^ 
As we will see later, this simplified treatment yields rea- 
sonable results in the present case as well. 

The staggered order parameter (^'^"^ (static value of 
yj) is zero below and above the critical field it is 
determined as the nontrivial solution of the equations 



(10) 



where the matrices Q, A are defined as 

hkhi 



kin 

A = r ■ -k r 



(11) 



lik 



kls 



' i? / , ^ikm^jln — — 

mirrij 



kl 



Linearizing the theory around ip = 93^''^ , one finds the 
magnon energies E depending on the field H and the 
wave vector q as real roots of the secular equation 

det(M -i;^G-i£;C) =0, (12) 

where = Qij{(fi^^'') and the matrix M is given by 



n. 



1 



^3 ' ■ ekY 

k 



(13) 



The antisymmetric matrix C can be written as 

c, = (i-+^)E'.,A-Eri"vr' (") 

'■ ■' k kl 

X {Sikij + Skii.j + Skiij — Sjki^i — Skji^i — Skij^i), 
where Sikij is defined as follows: 

SikLj = IkjSkl^^ e.ijrhr + ~\~'' f-ilrhr- (15) 
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To proceed, one has to fix the principal anisotropy 
axes. As suggested in Ref. ITol for symmetry reasons one 
of them should coincide with the crystallographic b axis; 
this axis we will denote z. Another one, which we label y, 
should be the axis along which the spin ordering occurs 
at Jf II b; according to Ref. 0, it lies in the (ac) plane 
and forms the angle of 39° with the a axis, see Fig. Q] 
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FIG. 2: ESR frequency-field dependencies!^ taken at T = 
1.5 K for different field orientations, in comparison with the 
present theory. Solid and dashed lines correspond to the 
ground state and thermally activated transitions, respectively. 



The best fit to the entire set of the experimental 
frequency-field dependencies is obtained with 

/3 = {5.628, 5.638, 5.635}, m = {5.725, 5.711, 5.714}. 

(16) 

Here all values are given in meV, the absolute error for 
the anisotropic part of /3, m is about 10~^, but for the 
isotropic part it is larger, about 0.1: the results are not 
very sensitive to a shift of (3i and fhi by the same value. 

We assumed the g-factor to be diagonal in the chosen 
anisotropy axes and took g^x — 2.29, gyy — gzz — 2.06. 
The theoretical curves are shown in Fig. [3 in compar- 
ison to the experimental results. All gaps occur at 
q = 0. Solid lines show the lowest magnon gap, and the 
dashed ones correspond to the transitions between three 
branches of the magnon triplet at g = 0. The theory is in 
a good agreement with the low-temperature ESR results. 
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FIG. 3: INS results of Rtiegg et alsM taken at T = 1.5 K, 
compared with the present theory. In this geometry (H is 
perpendicular to the plane defined by (010) and (104) vectors) 
the predicted lowest gap stays below 0.09 meV at H > He, 
i.e., beyond the available experimental resolution. 

We were not able to estimate ot from our fit, since the 
results are insensitive to the exact value of af. indeed, 
a enters only the 77^93^-type part of the interaction ||SJ), 
and in the ordered phase rj oc {h/J)ip, which suppresses 
the contribution of such terms in the relevant field range 
h/J< 0.1. Available estimatesi^ of J(ioo)i i) sug- 
gest that ai— pi « —0.92 meV, which has been assumed 
for the curves presented in Fig. |3 

Above the critical field the lowest gap opens as Eg ~ 
C\/h^ — /i^, e.g., for || 2; and weak anisotropy, 

~ 2{rhy/my){mymx - mxmy)/{rhx + rfiy)'^ ~ 0.16, 
where the weak contribution of 7^ , Ay is neglected. 



IV. DISCUSSION 

Our analysis suggests that the anisotropy in intradimer 
interactions rhi as well as in inter-dimer couplings [3i is 
very small and does not exceed one percent, which is 
plausible for the exchange anisotropy>i& From Eq. Hlfci|) 
one may get the impression that the intradimer exchange 
has a different sign of anisotropy. However, m given 
in (|16|l corresponds to the physical intradimer coupling 
J = {5.700, 5.728, 5.722}, so that the anisotropy has the 
same character for inter- and intradimer exchange: y is 
the easy axis, and z = h is the intermediate axis. 

According to our calculations, the staggered order pa- 
rameter <p(°) is directed along the y axis for H \\ b and 
H ± (102) (with a tiny deviation from the {yz) plane in 
the latter case); for H \\ [201], we predict that c^'"' || b. 
The magnitude of the order parameter at 7? = 90 kOe 
is about 0.1 of the saturation value, for all three field 
geometries. Those predictions agree with the existing re- 
sults of elastic neutron experiments^ for || 6, and can 
be tested for other orientations. Our results are consis- 
tent with the INS d&tdi^ as well, see Fig. 13 

It is worthwhile to note that the same conclusion on the 
character of the anisotropy {y and h being the easy and 
the second easy axis) was reached in recent ESR studiesi^ 
of the impurity- induced order in Mg-doped TlCuCla. 
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